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The possible role of the van Hove singularity (vHs) in stabilizing the low-temperature orthorhombic 
(LTO) phase transition in La2_a;Sra;Cu04 (LSCO) is discussed. It is found that the vHs can drive a structural 
' distortion in two different ways, either due to spin-orbit coupling or to dynamic Jahn- Teller (JT) effects. 

This paper discusses the latter effect in some detail. It is shown that a model Hamiltonian introduced earlier 
to describe the coupled electron - octahedral tilt motions ('cageons') has a series of phase transitions, from 
\ a high-temperature disordered JT phase (similar to the high-temperature tetragonal phase of LSCO) to an 

' intermediate temperature dynamic JT phase, of average orthorhombic symmetry (the LTO phase) to a low 

\ temperature static JT phase (the low temperature tetragonal phase). For some parameter values, the static 

' JT phase is absent. 

o 

cn 



o 



1. Introduction 



' The high-Tc superconductors La2_KBaa;Cu04 (LBCO) and La2-a;Sra;Cu04 (LSCO) have structural 

^ , phase transitions from a high-temperature tetragonal (HTT) to a low-temperature orthorhombic (LTO) to 

(in doped LBCO) a low-temperature tetragonal (LTT) phase[l] . There is a clear association of the LTT phase 
with electronic properties - the transition causes a decrease in the Hall density and interferes destructively 
with superconductivity [2]. Moreover, the transition is complete only at a fixed hole density, x 0.125[3]. 
Q . These features can be understood [4] in terms of the LTT phase splitting the degeneracy of the two van Hove 

singularities (vHs's)[5-7], at the X- and y-points of the (HTT) Brillouin zone. 

The role of electron-phonon coupling in the LTO phase is less clear, since transport properties are only 
' weakly affected by the HTT-LTO transition. While a uniform LTO phase could be stabilized by purely non- 

electronic mechanisms, this does not explain why high-Tc superconductivity appears to occur only in the 
orthorhombic phase[8]. It had early been suggested that the LTO phase was a charge-density wave (CDW) 
phase, associated with the vHs, but Pouget, et al.[9] suggested that this interpretation was untenable, since 
the LTO transition does not split the degeneracy of the two vHs's. In the present paper, I show that the 
analysis of Pouget, et al., is incomplete: the vHs can actually drive two transitions of macroscopic LTO 
symmetry, to either a CDW-like phase or to a spin-density wave (SDW)-like phase. 

I have recently suggested[10,ll] that both LTT and LTO phases are manifestations of a novel form of 
band Jahn- Teller (JT) effect[12-14], in which the degenerate electronic states are associated with the two 
vHs's. Thus, the LTT phase involves an essentially static band JT distortion, splitting the degeneracy of 
the two vHs's. The HTT and LTO phases can then both be interpreted as dynamic Jahn- Teller phases, 
involving tunneling between the X- and Y- point JT distortions of the LTT phase. In Ref. [10], a mean 
field calculation was made, approximating the dynamic JT effect by 'valence bond density waves', a coherent 
superposition of two charge density waves (CDW's). The resulting phase diagram reproduced the HTT 
LTO superconducting phases as a function of hole doping in LSCO. Ref [11] introduced a model 
Hamiltonian to describe the (nonlinear) electron-phonon interaction, based on similar calculations for the 
A15 compounds [15, 16], and applied it to an analysis of the static JT effect. The present paper extends these 
calculations to include the dynamic JT effect. For a single cell, the combined electron-tilt phonon ('cageon') 
problem can be reduced to the problem of a particle moving among four potential wells, and a solution can be 
found in terms of Mathieu's functions. The ground states are linear combinations of tilted octahedra. Adding 
intercell coupling leads to phase transitions with a net macroscopic average tilt, and the accompanying static 
strains. A dynamic LTO phase can be stabilized over a considerable temperature range by entropic effects, 
and the sequence HTT^LTO^LTT is naturally reproduced. A sequence of transitions is found, from a high 
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temperature disordered JT phase (equivalent to the HTT phase of LSCO) to an intermediate dynamic JT 
phase, with equal tilts along the orthogonal x and y axes (an LTO phase) to a low temperature static JT 
phase (LTT phase). The transition to the LTO phase is second order, while that between LTO and LTT is 
first order. Within a certain parameter range, the LTT phase does not occur, leaving the dynamic JT phase 
as the stable low-T phase. 

The dynamic JT phase offers an interpretation of the LTO phase as having only macroscopic average 
orthorhombic symmetry, with local dynamic disorder. The question of whether the vHs can induce a tran- 
sition into a uniform LTO phase is also reanalyzed, on the basis of group theory. It is found that umklapp 
scattering can split the vHs degeneracy, but only in the presence of spin-orbit coupling. Thus there are two 
vHs-based mechanisms of driving the LTO transition, in competition with each other. While considerably 
more complicated, these transitions are reminiscent of the competition between CDW and SDW transitions 
in the theory of nesting instabilities in lower- dimensional metals. 

The paper is organized as follows. Section 2 reviews the proposed role of band JT effects in the cuprates, 
and their relation to ferroelectric perovskites. Section 3 introduces a modified form of the electron-tilt-strain 
coupled Hamiltonian of Ref. [11], including umklapp scattering terms. Section 4 and Appendix II show 
how the acoustic strains can be formally eliminated, leading to an effective electron-electron interaction. 
Section 5 and Appendix III carefully analyze the problem of structural phase transitions in the uniform LTO 
phase. It is shown that, whereas umklapp scattering can in principle lead to a density wave instability, the 
transition is symmetry forbidden for the LTO space group, imless spin-orbit interaction is included. Section 
6 analyzes the dynamic JT state, reducing the Hamiltonian to an intracell plus near neighbor coupling form, 
and shows the relation between the present Hamiltonian and previous JT calculations. Section 7 presents the 
calculation of the dynamic JT effects at a mean-field level, including phase diagrams of the transitions from 
HTT^LTO^LTT. A discussion is given in Section 8, while the interpretation of the 'cageon' in terms of 
polarons or solitons is briefly discussed in Section 9. A number of Appendices provide details of calculations: 
a renormalization of the band structure to reveal separate subbands associated with each vHs (Appendix 
I); a group theoretical analysis of the uniform LTO phase (Appendix HI); a pseudospin approximation of 
the electronic operators (Appendix IV); and a discussion of the solutions of the four- well Mathieu equation 
(Appendix V). 

2. First-Order Band JT Effect 

The JT theorem states that any orbital electronic degeneracy in a molecule is unstable: there always 
exists some structural distortion which lowers the energy of the molecule by lifting the degeneracy. A similar 
effect can arise in a solid, either for a local impurity or for a collective JT effect in the crystal as a whole. 
Many perovskites and related structures have a phase transition to either a ferroelectric or antiferrodistortive 
phase, which can be interpreted as a collective JT effect[13]. 

In the cuprates, the structural transitions are complicated, involving tilts of the oxygen octahedra 
(optical phonons), static strains and acoustic phonons, as well as electron-phonon coupling. Moreover, 
the coupling of the optical phonons to both the strains and the electrons is nonlinear. A model for these 
interactions was introduced in Ref. [11]. There, it was suggested that both the LTO and LTT phase 
transitions were driven by the diverging electronic susceptibility associated with the vHs. Indeed, these 
transitions can be considered as a generalization to two dimensitions of the Peierls instability. 

The LTT phase transition can be thought of as a vHs-JT transition. The presence of two vHs's, at the X 
and Y points of the Brillouin zone, provides the electronic degeneracy, while the tilting of the octahedra split 
this degeneracy, as in a traditional (band) JT effect. Interpretation of the LTO phase is more complicated, 
since the associated octahedral tilting leaves the vHs's degenerate, essentially because the two in-plane Cu-0 
distances are equal in the LTO phase[9]. However, this rules out only the simplest model of the CDW/SDW 
phase. In the present paper, I show that there remain a large number of potential ground states, stabilized 
by (partial) splitting of the vHs degeneracy. I suggest that two of these states may be actually realized in 
LSCO: a static, SDW-like phase near half filling (Section 5), and a dynamic, CDW-like phase in the doped 
material (Section 7). 

The difference between the two transitions can be understood from Fig. 1, which illustrates two different 
means of splitting the vHs peak in the dos. This Figure illustrates the Fermi surfaces corresponding to the 
LTT phase (Fig. la) and to a possible LTO phase (Fig. lb). The Fermi surfaces at the vHs have been 
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distorted from square^ by inchiding a finite oxygon-oxygen liopping energy, too in the dispersion relations. 

Figure la shows Fermi surfaces of the LTT phase for two different dopings - the Fermi level coincides 
with the X-point vHs at one doping, and with the F-point vHs at the other. In the LTT phase, the two 
Cu-0 distances arc no longer degenerate, so the Cu-0 hopping parameter becomes anisotropic, splitting the 
vHs degeneracy [4]. (Due to the unusual layering, this compound has a complicated c-axis dispersion: the 
roles of X-point and F-point vHs's are interchanged on every other layer. It might be thought that such 
a large c-axis dispersion would wash out the effects of the vHs, but in fact, it is a consequence of the vHs 
splitting in the LTT phase.) In the dynamic model of the LTO phase, the symmetry is microscopically LTT, 
thereby locally splitting the vHs degeneracy, but with a macroscopic orthorhombic strain. 

This is to be contrasted with the situation in the miiform LTO phase, Fig. lb. The LTO transition 
doubles the unit cell volume, thereby reducing the Brillouin zone volume by the same factor, and introducing 
a two-fold degeneracy in all bands. Unklapp scattering can in principle lift the degeneracy of the two bands, 
so that their respective vHs's coincide with the Fermi level at two different dopings, as illustrated in Fig. lb. 
Note that in Fig. la, the Fermi surface of a single band is illustrated at two different dopings, whereas Fig. 
lb shows the Fermi surfaces of two different bands at a single doping level. Note further that, within either 
band in Fig. lb, the two vHs's remain degenerate [9]. The splitting of Fig. lb can be used, for instance, 
to describe inter layer coupling in the cuprates. In this case, the two 'bands' would correspond to the Fermi 
surfaces at kz = 0,7r/c, which are degenerate when interlayer hopping is neglected. However, interlayer 
coupling is weak in these materials, and for simplicity it will be neglected. 

Of more interest is the possibility of introducing a gap by coherently coupling excitations between the 
two vHs's. Such in-plane coupling was introduced by Schulz[17], and is very similar to that introduced in the 
flux phase[18] of a doped antiferromagnct. However, a group theoretical analysis rules out such a possibility 
for a CDW-like excitation: there is a two-fold degeneracy of electronic levels on the Brillouin zone face, 
which can only be split by spin-orbit coupling. This suggests that vHs splitting can lead to a uniform LTO 
phase only in the presence of antiferromagnetic order - i.e., close to half filling. Away from half filling, the 
magnetic state becomes inhomogeneous[19]. Hence, in the doping range associated with superconductivity, 
there is competition between a disordered magnetic state and a dynamic JT state. Both phases will be 
described in greater detail below. 

It seems likely that dynamic effects are important in the higher-Tg phases, and that the sequence 
HTT— >LTO— >LTT may parallel the sequence of structural transitions in the perovskite, BaTiOs [20] . Here, 
four successive phase transitions are observed as temperature is lowered, from cubic to tetragonal to or- 
thorhombic to rhombohedral. The rhombohedral phase can be identified as a static JT phase, in which 
the Ti are all distorted along one of the eight octahedral directions. The higher T phases correspond to 
dynamic JT phases, in which the Ti tunnels between successively two, four, or all eight octahedral sites. 
In particular, the 'cubic' phase is not microscopically of cubic symmetry, but is a disordered phase with 
microscopic rhombohedral and macroscopic average cubic symmetry. A similar model can apply to LSCO, 
with a disordered or dynamic JT phase consisting of a random tilting of the oxygen octahedra about each 
Cu, with tunneling between all four allowed tilts in the HTT phase, between two in the LTO phase, and 
with only one tilt (static JT) in the LTT phase. Section 7 shows that such a sequence of phase transitions 
naturally follows from the Hamiltonian introduced[ll] to describe this system. 

Recent calculations[ll,21] have confirmed the close relation between the ferroelectric perovskites and 
the superconducting cuprates. In [11] it was pointed out that electron-phonon interaction can lead to a 
structural instability (negative harmonic phonon frequency) even in a filled band, as long as the Fermi level 
fell between the bonding and antibonding bands associated with hybridizing atoms. This is clearly the case 
in La2Cu04, and Cohen[21] has shown that Ti-O hybridization is essential for ferroelectricity in BaTiOs and 
PbTiOs. 

3. Theory of the Structural Transition 

The present calculations are based on the nonlinear clcctron-phonon Hamiltonian introduced in Ref. 
[11]. For convenience, when referring to equations from References [10] or [11], I will prefix them with a 
letter A ov B respectively (i.e., Eq. 10 of Ref. [11] will be referred to as Eq. BIO, etc.). 

It is convenient to rewrite the Hamiltonian of Ref. [11] as 

H = Hstr + Hph + Hf, + Hps + Hes + Hep. (1) 
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The various terms of Eq. 1 refer to 

• static strains and acoustic phonons 

= ^ + Y.{c+e+{-q)e+{q) + C_e_(-g)e_(9) + 2Ceeei2{-q)eu{q)) ■ (2) 

Xq q 

Here, Cij is the usual strain tensor, the elastic constants are Cij[6], and e± = (en ± e22)/2, C± = Cn ± €22- 
The strain tensor includes both static and dynamic parts (both strains and acoustic phonons), and the 
P(A,(j')'s are the momenta associated with the acoustic phonons. In the HTT phase, there may be static 
strains due to (anisotropic) thermal expansion: en ^ 0. A non- vanishing e\2 arises in the LTO phase. In the 
HTT and LTO phases, en = 622- This is also true in the LTT phase of LBCO, due to intcrlaycr strains[ll], 
but need not be true in general. To simplify the present calculations, I will apply them to a 'single layer' 
model of LBCO, for which e_ ^ in the LTT phase. 

• electronic modes 

In the vHs model, there is a single, hybridized Cu-0 band at the Fermi level. However, as discussed 
in Section 2, it is convenient to treat the two vHs as independent, since they have opposite responses to 
e_ strains. Hence, the carriers will be assumed to belong to group 1 or 2, depending on whether they are 
nearer the X- or Y-point vHs, respectively. In Appendix I, this assumption will be justified by calculating 
such subbands from a renormalized microscopic Hamiltonian of the hybridized Cu-0 planes. The electronic 
energy may be written as 

H, = J2Mk)4k<^ik, (3) 

i,k 

with i=l,2. It is convenient to define the integrated density of states for each vHs as 

Pijil) = ^<4kO'j,k+q>, (4a) 

k 

P± (q) = Pii (q) ± P22 {q) ■ (46) 

In the HTT and LTO phases, pu = P22, whereas they are unequal (the JT splitting) in the LTT phase. 
The term involving pi2, which is in general allowed by symmetry, corresponds to inter- vHs coupling. This 
term is related to umklapp scattering between the two vHs, and can act as a competing mechanism to the 
dynamic JT effect. In Section 5 it will be shown that such a term could drive a structural transition within 
a purely orthorhombic phase. However, it will further be shown that pu is symmetry forbidden at the vHs 
in the uniform LTO phase, in the absence of spin-orbit coupling. 

In the present paper, only intraband inter-vHs coupling will be considered, with pi2 coupling the vHs's 
separated by wave vector Qo = {n/a,TT/a), as in Refs. [10,11]. It should, however, be noted that similar 
considerations could be applied to interlayer coupling, as discussed in Section 2 (Fig. lb). 

• optical (tilting) modes 

Tilting of the oxygen octahedra can be described in terms of (pseudo) rotation operators, for rotations 

about the in-planc x- and y-axes[ll,22], with i?,. and Ry being the magnitude of the tilt. To simplify the 
resulting expressions, it is convenient to define some auxilliary pair tilting operators. Thus, 

Rij{q, q') =Ri{q + q')Rj{q - q'), (5a) 

nf^{q) = Y.RiA<l^q'), (56) 
q' 

T^ifil) = ^RiM,<l')^i3{cos{qja) - cos{q'ja)), (5c) 
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with Sij =1 — 6ij = 1 when i ^ j, i,j = 1, 2. In this case, Hph becomes[ll] 

^J"^ = ^E(E(^mm(9)^/^(-9)^'m(9)) +'^o(9)(^ii(0,a) + ^^22(0,?)) 

q n 

+ y [^li^ (3)^1 (-3) + (-«)] 

where P/j, is the momentum conjugate to i?^, /x runs over 1 and 2 (a; and y), and i^T^^u is defined in Eq. B17a. 

• electron-strain coupling The strain-electronic interaction is 

Hes = Y,{G+e+{-q)p+{q) + G.e.{-q)p.{q) + 2G^eei2{-q)pi2{q)) ■ (6) 

with G± — G\i± G22- The term in Gee couples the orthorhombic shear strain ei2 of the LTO phase to the 
umklapp term, < pi2 >. 

• strain-tilt coupling 

The optical phonon coupling with strain (Hps) can be written 

Hp, = ^(F+e+(-g)7ef (|) +F_e_(-9)7eL')(|) +F66ei2(-g)7e^'^(|)). (7a) 

Here, I have introduced the combination 

Tef {q) = cos{qya)TZf^ {q) ± cos{q,a)TZ^^ {q) , (76) 

and 

T^o\l) = ^Ri2{q,q'){l - cos{{qy + q'y)a))il - cos{q, - q'Ja)). (7c) 

9' 

These equations reduce to the corresponding expressions in Ref. [11] when q = 0. In that reference, the 
Hamiltonian was evaluated only for wave numbers near the soft mode, q = 0, q' at one of the vHs (either 
(0,7r/a) or (7r/a,0)). 

• electron-tilt coupling The optical phonon coupling with electrons {H^p) can be written 

H,p = ^(a;p+(-«)7ef (|) + a^_p_(-g)7eL')(|) +7Vi2(-9)7e^'^(|)). (8) 

Q 

In the notation of Ref. [11], aj_ = {S'^ ± a'^)/2m, 7^ = j^/2m. This is a nonlinear electron-phonon coupling, 
similar to terms which have recently been introduced [23]. 

All of the above terms involve a sum over q. In an RPA approximation, the various g's are decoupled. 
The soft mode is associated with g = 0, so in many applications the non-0 g'-terms may be neglected (e.g., 
only terms in 7?.j_ (0) remain). Moreover, except for the tilt term, Eq. 5d, the Hamiltonian separates into 

three parts, involving (e+, p+, TZ'f^), (e_, p-, 7?.f^^), and (ei2, P12, Ti-'o^)- In Appendix II, a toy model is 
introduced which allows a similar separation of the tilt terms, Eq. 5d, as well. 

The three groups of terms play very different roles in the structural transitions, particularly in the 
presence of dynamic JT effects. At high temperatures, corresponding to the HTT phase, there are local tilts of 
the octahedra, but with no long-range correlations. In this case, < R'^+Hy >^ 0, while < Rx >=< Ry >= 0. 
Hence, the first set of terms controls short range order, and is non- vanishing in the HTT phase. These terms 
will have a weak temperature dependence at lower temperatures, which can be neglected to simplify the 
study of the dynamics - i.e., the first group of terms contributes a constant value to the Hamiltonian and 
can be neglected. 
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The second group of terms is nonzero in the LTT phase, and the last in the LTO phase. Hence, these 
terms describe the dynamic competition between the LTT and LTO phases. 



4. Eliminating the Strain and Acoustic Phonon Modes 

The elastic strain tensor has static components associated with static strains and time-dependent 
components associated with acoustic phonons[ll,22]. However, in studying structural phase transitions, it 
is often convenient to treat the strains separately from the acoustic phonons. For instance, condensation of 
an optical mode may induce a static strain. Consider a line of corner-shared octahedra: if the interatomic 
distances remain fixed, a static, antiferrodistortive tilt distortion of the octahedra will reduce the overall 
length of the chain. This seems to be the case in the cupratcs: the octahedral strain in the LTO phase is 
found to be a secondary order parameter, proportional to the square of the octahedral tilt angle [24]. These 
strains may be formally decoupled from the problem, following standard practice[12-14], by defining 

6p — 6-p ~\~ Ap, (9) 

(p = —,12) and chosing Ap in such a way as to eliminate the cross terms between Cp and pp, RiRj. The 
procedure is carried out in detail in Appendix H, but for illustrative purposes, a simpler calculation is given 
here. Consider the sub-Hamiltonian 

H' = 2C66e?2 + 2G66ei2Pi2. (10a) 
Taking A12 = G66P12/2C66 transforms the Hamiltonian to 

H' = 2C66gL - 4^P?2, (106) 

with J12 = Gqq/Cqq. Thus, the phonon motion is formally decoupled from the electron, leading to an 
attractive interaction between electrons. However, an electronic phase transition is accompanied by a static 
distortion: since when < 612 >= 0, then 

< ei2 >= < P12 > . (11) 

The separation of charge and phonon variables is not complete (ei2 does not commute with pu), but it 
has been argued that the additional complications of noncommutivity arc unimportant (see discussion in 
Ref. [14], pp. 24-25, and references cited therein). The resulting attractive interaction between electrons is 
very similar to that found in the density wave calculation [10, 2 5], Eq. 2.8 of Ref. [25], which, in the present 
notation, becomes 

Y ^ GlQhujq/2 GIq 

Equations 11 and 12 difii'er only in the denominator, with one equation containing Cee, the other Uuiq. This 
substitution arises quite naturally. In Eq. 2, the term 612(5) includes both static strains and acoustic 
phonons. If the phonons only were included, the term 2C66ei2(— (z)ei2(?) cl^Cqihojq + 1/2), where the Cq 
are phonon operators. Elimination of the electron-phonon coupling as above would then lead to an effective 
electron-electron coupling, Eq. 11 with Cge fiujq. 

In Appendix II, a toy model Hamiltonian is introduced, which allows the optical phonon coupling to be 
eliminated in the same fashion, leaving a purely electronic Hamiltonian with an attractive effective electron- 
electron interaction. There is also the decoupled phonon Hamiltonian, which splits into two parts: (1) an 
acoustic phonon part which, being purely harmonic, can be neglected; (2) an optical phonon part, which 
remains anharmonic, but with renormalized coefficients. 

5. Uniform LTO Phase 

Before analyzing the dynamic JT phase, it is important to reexamine the question of whether a structural 
transition from the HTT phase to a uniform phase of microscopic LTO symmetry could somehow be driven 
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by the vHs. At first sight, it would appear that this question has already been answered in the negative by 
the work of Pouget, et al.[9]. However, their analysis did not consider all possible mechanisms for driving a 
structural instability. 

Elcctron-phonon interaction can lead to a structural instability if the structural distortion drives a 
significant density of electronic states below the Fermi level, as in the one- dimensional CDW problem, due 
to Fermi surface nesting. In the presence of two degenerate vHs, there are two different ways in which a large 
dos could be shifted below the Fermi level. First, the structural distortion could split the degeneracy of the 
two vHs, so that one vHs is shifted below the Fermi level, the other above. This is clearly what happens in 
the LTT phase of LBCO, Fig. la, and is the basis for the dynamic JT model of the LTO phase, discussed in 
Sections 6 and 7. As shown by Pouget, et al.[9], the vHs degeneracy is not lifted in the uniform LTO phase. 

However, there is an alternate mechanism for structural distortion, which could in principle provide a 
static model for a uniform LTO phase. In this mechanism, umklapp scattering couples electronic states on 
degenerate bands of the Fermi surface. Fig. lb, leading to coherent superpositions of the two states with 
a corresponding gap between the superposed states. In the LTO phase, there is no splitting of the vHs 
(Fig. lb), but the two vHs are now at symmetry-equivalent points of the Brillouin zone, thereby allowing 
inter-vHs umklapp scattering. By introducing a gap at the vHs, this mechanism can again stabilize a lattice 
instability. In the Equations, 1-8, this inter-vHs coherence corresponds to a finite expectation value for pi2, 
and the umklapp scattering to the Gee and 7^ terms. For the calculations of this section, it will be assumed 
that the two vHs's are separated by Qo = (7r/a,7r/a), although, as discussed below Eq. 4, other Q-vectors 
are possible. 

The present section is arranged as follows. In Section 5a, the umklapp mechanism is discussed, and it 

is shown that it is likely to be weak in doped LSCO: it is symmetry- forbidden in the absence of spin-orbit 
interaction. Section 5b will further demonstrate that the static model cannot explain the sequence of transi- 
tions from HTT to LTO to LTT in terms of competition between pi2 and p- . There is indeed competition, 
but the present calculations suggest an either/or situation: if there is a static structural instability, it will 
be either to an LTO phase or to the LTT phase, depending on the relative parameter values. Thus, the 
experimentally observed sequence of phases suggests rather a dynamic LTO phase - as in the perovskite 
ferroelectrics. The calculations of Sections 6 and 7 confirm this possibility. The correct sequence of phase 
transitions is found, even though umklapp scattering is explicitly neglected (Gee and 7'^ set equal to zero). 

5a. Interpretation of pi2 

The significance of the term pi2, Eq. 4, can be clarified by recalling the usual manner in which gaps in 
the electronic spectrum open at a Brillouin zone boundary [26]. Umklapp scattering mixes states at k and 
k + Q. Thus, the electron at the zone boundary k = Q/2 mixes with that at A; = —Qjl, 

il) = a+cl^^2 + a_cLQ/2' (13a) 
with the a's determined by the eigenvalue equation 



{Eq/2 - E)a+ + Uga- = 0, (136) 

UQa+ + {E_Q,2 - E)a^ = 0, (13c) 

where Uq is a measure of the coupling. Now in the LTO phase, the X and Y point vHs's are separated 
by a reciprocal lattice vector, so a similar coupling can arise. In the present formalism, this coupling arises 
by rederiving Eq. B24 in the presence of the term in 612^12- Neglecting fluctuating quantities, Eq. B24 

becomes 

d ~ ~* - 

i^ank = En{k)ank + E*amkSmn, (14) 
ot 

Ei{k)=Ei{k) + G-e- + a'LRl, (15a) 
E2{k) = E2{k) - G-e- - atR^_, (156) 
E*{k) = 2Geeei2 + ^^RiR2. (15c) 
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Equation 14 shows that, because of the /9i2-term, carriers from the two vHs interact with each other, in 
the presence of an orthorhombic strain (ei2 ^ or i?i = ±i?2 7^ 0). In this case, the correct eigenstates can 
be found by the BogoHubov construction. Defining 

aik = cos{0)cik + sin{9)c2k, (16a) 

a2k = -sin{6)cik + cos{6)c2k, (166) 
then the equations of motion for the c,fe's are decoupled if 

E2{k)-Ei{k) 

and the eigenenergies become 



tan{2e) = - , (17a) 



E*^=E+± y/sl+E*'^, (176) 

with E± = {El ± E2)/2. The term E_^. may be neglected in Eq. 17b, since any term which shifts both 
electronic bands equally will be compensated by a corresponding shift of the Fermi level. Thus, when a 
static orthorhombic strain appears (< ei2 >^ 0), a gap 2E^ arises in the electron spectrum, driving the 
high density of states associated with the vHs below the Fermi level. 

In fact, however, the terms in p-12 must vanish identically in a uniform LTO phase, in the absence of 
spin-orbit interaction. This follows from the symmetry group Bmab of the LTO phase, and can be considered 
as a generalization of Pouget, et al.'s result[9]. Because the lattice contains glide planes, the eigenfunctions 
must be two-fold degenerate on one face of the Brillouin zone, so no gap (pi2 ^ 0) can arise. However, 
spin-orbit interaction splits most of the residual degeneracy, allowing a gap to open at the vHs. The group 
theoretical arguments are discussed in more detail in Appendix III, and the effect of such spin-orbit coupling 
on the vHs will be discussed in Section 8. 

5b. Mean Field Transition in Electronic Hamiltonian 

The above formalism can also be used to study the competition between static distortions of LTO vs 
LTT symmetry. This is most clearly seen by analyzing the effective electron-electron coupling terms in Eq. 
10b. To simplify this analysis, it is convenient to temporarily neglect the tilt coupling and study just the 
electronic Hamiltonian. Alternately, the toy Hamiltonian of Appendix II can be used to formally eliminate 
the tilt-electron coupling. From Appendix II the effective electronic Hamiltonian is found to be 
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i ^(j_p_(-q)p_(g) + JoPi2(-g)pi2(<z)) , (18) 



Q 



with the coupling constants, Ji, defined in Eq. 116. A term in p+ has been neglected in Eq. 18. From charge 
conservation, p+(0) must be a fixed constant, which can be set equal to zero by adjusting the Fermi level. 
Due to the logarithmic divergence of the electronic susceptibility at the vHs, the dominant singularity (soft 
mode) corresponds to g = 0, and at the RPA level of approximation, this is the only mode which need be 
discussed. 

The mean field solution can be found easily, as in Section 5a and Ref. [10]. The mean field Hamiltonian 
becomes 

Hmf = X^(-£^fc(«ifcaifc + a\f.a2k) - D_{a\^aik - oL«2fc) - Do{a\^a2k + alk<^ik)^, (19a) 
k 

with 

£>_ = ^ <p_(0)>, (196) 
^0 = Y < Pi2(0) > . (19c) 
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H%fp can be diagonalized as in Eqs. 16, 17, yielding eigenvalues 

El^ =Ek±D, (20) 
with Z)^ = D"^ + Dq. The self-consistency conditions, Eqs. 19b, c, yield equations for the two gaps 

fe 

fe 

where f{E) = l/{exp{{E — Ep)/kBT) + 1). Except in the special case J_ = Jq, Eqs. 21a, b cannot 
simultaneously have nonvanishing solutions. Instead, the solution with the larger value of Jj prevails. This 
makes good physical sense: the two transitions, LTO and LTT, are both driven by the same dos peak. 
Whichever phase is stabilized first uses up the available dos, and prevents the other from occuring. 

In this case, the calculation may be simplified by ignoring terms of the non- condensing symmetry. This 
was in effect what was done in Rcfs. [10] and [11]. In Ref. [10], only the LTO solution appeared, while 
the calculations of Ref. [11] neglected terms involving pi2, and predominantly described the LTT phase. It 
should be particularly noted that the present calculation has reproduced the BCS-like calculations of Ref. 
[10], but starting from the more microscopic Hamiltonian of Ref. [11], thereby explicitly displaying the close 
connection between the two works. 

To compare with the results of Ref. [10], assume the LTO phase is favored, i.e., Jq > J_. Then £)_ = 0, 
and the gap D = Dq is given by the solution of 

1 = ^ / dENiE) {f{E-D)-f{E + D)). (22) 

For a logarithmic dos, 

N(FA = 

B ^2E 

the zero temperature gap is 



N{E) = hn{^), (23) 



i)(0) = ^e-^/^«. (24) 

This should be compared to Eq. 14 of Ref. [26]. Figure 2 illustrates the temperature dependence of the 
gap, found by solving Eq. 22-23 numerically. Note that, based on Eq. 11, there will be a nonvanishing 
orthorhombic strain in the LTO phase, < ei2 >oc — < pi2 >. Indeed, the present solution is similar to that 
found in Ref. [11], Model 2, except that inclusion of umklapp processes drives the large vHs dos below the 
Fermi level. 

This section has explored the role of the term pi2 in stabilizing a static LTO phase. However, symmetry 
arguments suggest that such terms arc small in the absence of spin-orbit coupling. Hence, in Sections 6 
and 7, the opposite limit will be explored. Terms in pi2 will formally be retained in the Hamiltonian, but 
I will attempt to determine under what circumstances an LTO-type phase might arise when the pu terms 
are small or vanishing. It will be shown that the LTO phase can be interpreted as a dynamic JT phase. 



6. Dynamic JT Hamiltonian 



6a. Real-Space Hamiltonian: Intracell and Intercell Coupling 

This section will present a more accurate treatment of the combined electron-optical phonon Hamilto- 
nian, following a conventional treatment of the dynamic JT effect. It is convenient to first transform the 
Hamiltonian, Eq. 1, back into real space: 

l,ri 
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+{atp-{i) + i^_e_(r))[(Ml:-|te^)2 _ (MbMjlZ^)^] 

^^ ^^12(0+^66612(0 ^ ^[i?i(r) - R,{1 + wa)][i?2(Z) - R2{1 + v' xa)] 

v' 

+ '^{Ri{lf + R2{lf) + C+el{l) + C-el{l) + 2C66eL(0 

+G+e+{l)p+{f) + G_e_(r)p_(0 + 2G66ei2(0pi2(0) , (25) 
where rj and 77' are summed over ±1 and 

p-{l) = a\iau - a\ia2u (26) 

with a similar expression for p\2- To form the full Hamiltonian, Eq. 1, the electronic energy, Eq. 3, and the 

kinetic energy terms of Eqs. 2 and 5d must be added to Eq. 25. Moreover, a possible g-dependence of wq 
has been neglected. Eliminating the strain terms from Eq. 25, as in Appendix II, yields 

l,ri 

I a%p^{V)[{ ^^^^^ ~ ^ ^"'^'^ f + ^ R^jl) - R2il + yax) 

^ji^^^ Rijl) - RiQ + rjay) 2 , R2{l)-R2{l + vax) ^2^ 
+a_p-{l)[{ ) -( )] 

+^pi2{i) YliMi) - vya)][R2{f) - R2{1 + v'xa)] 



+^{Rr{lf + R2{lf) - ^-fplil) - ^p'.if) - f p?2(r)) + H*, (27a) 

where ^ = Eq - FI/C+ - Fl/C, a% =a%- G±F±/4C±, = 7' - 066^^66/2^66, J'± = Gl/2C±, and 
Jg = Ggg/Gee; H^tr is a quadratic pseudostrain Hamiltonian (Eq. Hal) decoupled from the remaining terms 
and of no further interest; and 

= -^E(ro+[(^i(0 - RiiU ay)f{R^{l) - R,{1- ay)f + 

{R2{1) ~ R2{1 + ax)f{R2{l) - R2{1- ax)f] 
+ro_[(i?i(r) - Ri{l + ay)f + {Ri{f) - Ri{l - ay)f][{R2{l) - R2{1 + ax)f + {R2{1) - R2{1- ax)f] 

+roo(53[i?i(r) - Ri{l + vya)][R2{l) - R2{1 + 'n'xa)]f), (276) 

with Toi = F2/G+±F2/G_, Too = ^16/8^66- 

The role of the strain forces can bo determined by comparing Eqs. 25 and 27. In addition to the 
effective electron-electron interaction terms, the strain has introduced longer-range tilt-tilt interactions. All 
of the terms of Eq. 25 are either on-site or nearest neighbor interactions, except the terms in 7* and Fee- 
In contrast, all of the terms of H* involve further neighbor interactions. 
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In a mean field treatment, it is assumed that there is a nonzero tilt present on each lattice site, even in 
the HTT phase 

< Rfil) + Rjil) >= i- 0, (28a) 
independent of /. Then the dynamic variable is the tilt direction, 0;, with 

iJi(r) = {-ly^^ Rcos^u (286) 

J?2(r) = (-Vf"^^Rsin^u (28c) 

with I = {ia,ja). Note the factor (—1)*+-' - this is introduced because the intercell coupling must locally be 
antiferrodistortive, due to the corner sharing of the octahedra. By explicitly taking out this factor, it can 
be expected that (pi will be a smooth function of position, with a well behaved continuum limit. 

Equation 27 can be separated into intracell and intercell terms, as H = Xl((^o + + ^2): ^it^i 



(29a) 



i - of' - 'Y*' 

= --T'^R'^cosAcj) + -^p_R^cos2(j) + ^ pi2R'^ sin2(f) 



-^pi - (296) 

with rff = rg + (r;,-ro+)/16 and r^' = (r2-2r;,'-roo-ro-/4)/8. since all terms refer to the same cell, the 
/-dependence is not explicitly displayed. Here, Hg is 0-independent, is a single cell Hamiltonian, while 
H2 incorporates the intercell coupling. The intercell coupling term is complicated, and will be explicitly 
displayed only at the mean field level, for which cos<j){l') =< coscf) > and sin(j){l') =< sincf) > for all I' 7^ I. 
In this case, iJf — -f^2(< coscj) >, < sin0 >) — i?2(0) 0)) with 

H2{< coscf) >, < sinct) >) = ° R'^[{cos(t)+ < coscj) >)'' + (szn(^-|- < sintt) >)^] 

— ^^^'^ R^(cos(j>+ < cos(f> >)'^{sin(p+ < sincf) >)^ 



H — — [{coscj)+ < coscj) >)^ + {sincj)+ < sincf) >)^] 



p R? 

H — ^—^ \(coscf)+ < coscf) >Y' — {sincf)+ < sincf) >)^] 

+^'" p\2R^{coscf)+ < coscf) >){sincf)+ < sincf) >). (29c) 

As discussed in [11], the terms in /?+ are non-critical, and can be eliminated from Eq. 27. Thus, 
p+{q = 0) is just the number of holes in the conduction band, p+(0) = 1 in the present case, where P-(0) 
(or P12) becomes non-zero only in the low-T phase, and hence can be taken as an order parameter of the 
transition. Thus, we may assume p+{q) ~ P+{0), and eliminate the p+-dependent terms from Eq. 27. The 
harmonic phonon frequency is renormalized u)q = ojQ + a^{p^{0) — 2). In this equation, I have incorporated 
an additional correction[ll], the term in —2, due to the filled, bonding band of the hybridized Cu-0 band. 
This term is important in destabilizing the lattice, u)q < 0, both in the present problem and in ferroelectrics 
and other structural instabilities. 

While the p+ -terms are neglected in the present analysis, they may yet have important effects in these 
materials. They provide a coupling between the tilt and the local average electronic density, and lead to the 
possibility of a microscopically heterogeneous phase. Such nanoscopic disorder has previously been suggested 
to play an important role in doping these materials away from the vHs[28] , while a related phase heterogeneity 
has been proposed to arise on doping away from the antiferromagnetic phase at half-filling[17]. 
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6b. Intracell Hamiltonian: Static JT Effect 

It is convenient to begin by discussing the JT effect within a single celL The terms in Hq lead to a 
tilting of the octahedron, i? 7^ 0, but with no preferred orientation. If R is assumed to be a fixed constant, 
then the JT effect involves the angular orientation of the tilted octahedron, (p. For the single cell problem, 
this involves H", Eq. 29b. The solution to Eq. 29b depends on whether the system is in the static or 
dynamic JT limit. In the dynamic limit, the tilt kinetic energy operator associated with (f> must be added 
to Eq. 29b; in the opposite limit, the tilts are static, and the kinetic energy can be neglected. In this static 
case, the electronic operators can be diagonalized by a transformation similar to Eq. 16. Alternatively, a 
pseudospin formalism[12-14,29] (Appendix IV) can be employed. In the two-dimensional subspace spanned 
by the electronic operators a\{l), a^il), the p operators can be represented by Pauli matrices 

P- = o-z, 

P12 = <yx] 

hence, the p^-terms in Eq. 29b reduce to constants, and can be eliminated. Transforming the electronic 
states by Eq. 16, the eigenenergies of Eq. 29b become 

= - ir^'^''cos4(/) ± i ^Ja''l^R'^cos'^2(j) + Y''^R^sin^2(j). (30) 

Equation 30 demonstrates the JT effect in the present system. If E°_ were minimized with respect to R, 
the electronic term would always give rise to a non- vanishing JT distortion (^ 7^ 0); incorporation of other 
terms from the original Hamiltonian of order ~ R'^ would not change this result. However, since the electron- 
phonon coupling is quadratic in R, there are also harmonic terms in the Hamiltonian of order i?^, and a JT 
splitting will arise only if the coupling coefficient (e.g., a"') is large enough (or the harmonic coefficient is 
negative see [11]). 

Equation 30 represents a potential with four degenerate energy minima. Depending on the parameter 
values, these minima may lie either along the ±a; and ±?/ axes, or at 45° to this. It is expected that 
a!L' > 7*"; if 7^' = 0, the minima would lie along the principal axes when —2T2R? — < 0, and at 45° 
when the inequality is reversed. In the static limit, the tilted octahedron will be located at one of the four 
equivalent minima. In the absence of intercell coupling, there can be no macroscopic phase transition: at 
low temperatures the octahedra will be randomly distributed among all four minima; as the temperature is 
raised, the octahedra can hop among the various minima. 

In what follows, it will be convenient to approximate the lower JT solution, of Eq. 30 by a simpler 

form 

El_ c^Ea- EbCos4(f). (31) 

For instance, if 7'^' = 0, 2Ei, ~ |al'|-R^/2 + T'^R^. For nonzero 7*', the angle dependence is a function of 
— 7*'. The terms which comprise E^ tend to appear as the difference between two quantities (e.g., Fq, 
Q,e/ _ jBiy rpj^jg occurs because for an isolated free molecule the 'JT distortions' amount to a pure rotation 
of the molecule; the distortion is a solid-state effect, due purely to crystalline anisotropy (this fact appears 
to have been first pointed out in Ref. [30]). The constant Ea can be absorbed into the angle- independent 
part of the Hamiltonian, H^. 

6c. Dynamic JT Effect and Relation to Conventional (E (g) e) JT Effect 

When dynamic effects are important, the ionic kinetic energy operator from Eq. 5d must be retained 
in the Hamiltonian, Eq. 29b. The exact eigenstates no longer can be written in Born-Oppenheimer form, 
but are of the form 

*„(i?, 0) = XnliR, 4>)^-{R, ^) + Xn2{R, <^)^+(i?, <^), 

where 'il)± is the electronic wavefunction corresponding to the energy E^^, Eq. 30, and the Xni are wave- 
functions of the nuclear motions. In the limit of strong JT coupling, the upper JT level can be neglected, 
and a Born-Oppenheimer wave function is approximately recovered: 

^n = f{R)cl>n{4')X±{RA), 
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where x± tiio lower energy electronic wave function and the nuclear wave func;tion has been separated into 
a radial part, f{R), assumed to be approximately constant, and an angular part, 0„, which satisfies 

{jrr2 + pcosmi))(f>n = Er,(l>n, (32a) 



where the first term in Eq. 32a is the angular contribution to the tilt kinetic energy, and 

This distorted octahedron, strongly and nonlincarly coupling the tilt and the electronic state, constitutes 
the polaron of the present problem. For convenience, this tilt polaron will be referred to as a 'cageon'. 

Equation 32a is a form of Matthieu's equation[31-33] which often arises in the dynamic JT problem. 
The vibronic potential here has four minima, along the positive and negative x- and y-axcs (for /? > 0), 
corresponding to the four possible static JT tilts of the CuOe octahedron. When P is very large, the lowest- 
energy state is fourfold degenerate, corresponding to these static distortions, with a weak tunneling between 
the states. For smaller barriers, this quartet of states breaks up into a pair of doublets. The solutions of Eq. 
32a are discussed further in Appendix V. 

In the closely related S ® e JT eff'ect[12-14], a similar Mathieu's equation arises, but with three-fold 
degeneracy (corresponding to elongation of the octahedron along an x, y, or z axis). This problem is often 
simplified by approximating the intercell coupling by a quadratic form. In this case, the weak tunneling 
limit reduces to a three-states Potts model[34,35]. For the present, four-fold degenerate model, this would 
correspond either to a pseudo-spin 3/2 system, or to a four-states Potts model. 

7. Intercell Coupling 

7a. High- Temperature Limit (Disordered JT Phase) 

In the intercell coupling term, < sirKf) > and < cos(p > are independent variables, constrained by (< 
sincj) >)^ + (< cos(f) >)^ < 1. Two classes of solution arc of particular interest: < cos4> >^ 0, < sirup >= 0, 
the LTT solution, and < sincj) >=< coscp >^ 0, the LTO solution. 

The analysis is most straightforward in the weak tunneling limit, when only the four lowest-lying levels 
of each octahedron need be considered. In the high-temperatmc limit, < sincj) >=< coscp >= 0, and only 
the single cell Hamiltonian, Eq. 29a,b, need be considered. If tunneling is completely absent, the four wave 
functions are each localized in one of the potential minima of cosAcf). Near the minimum, the potential is 
quadratic in (/), and the wave functions arc well approximated by harmonic oscillator wave functions. In 
particular, the ground state will be approximately a Gaussian wave function, 

^ = /|e-^^^ (33) 

with V = \/2^. Inclusion of overlap between Gaussians centered on difii'erent wells splits the degeneracy. 

However, just as in the 3-well model, it is important to recognize that the total wave function is a combined 
electron-phonon wave function. The electronic wave functions are only symmetric under a At: rotation; 
x{(j) + 2ti) = — x(0)- Thus, the nuclear wave function must also satisfy 

<^„(.^ + 27r) = -</.„(</.), (34) 

so that the total wave function ^„ has 27r symmetry. In the n=3 case, this changes the sign of the overlap, 
thereby reversing the order of the levels[12,31,36]. This sign change is now recognized to be an example of 
Berry's phase[37]. 

In a perturbation calculation of the four- well problem, this means diagonalizing an 8x8 matrix. However, 
since the wavefunctions 5-8 are just the negative of wave functions 1-4, this immediately reduces to a 4x4 
matrix 

Fii -E H12-SE -{H12 - SE) ' 

H\2 — SE Hii — E H12 — SE I — n 

H12 - SE Hii ~E H12-SE ' ^ ' ^ ^ 

-{H12-SE) H12-SE Hii-E 
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with S the wave function overlap, H12 the nearest-neighbor Hamiltonian overlap, and Hn the diagonal 
Hamiltonian. For the present problem, these matrix elements are 



Hn = av{l - "^e-^'n, (36a) 

Fi2 = 5aKl + ^(e-'/"'- y)), (366) 

S = e-'^"'/^ (36c) 
with a = /2R^. This overlap splits the four degenerate levels into two pairs of levels, with energies 

= H21±^. (36.) 

The wave functions associated with Ej^ are in the subspace spanned by the vectors (l/\/2, 1, l/\/2, 0) and 

(0, 1/V2, 1, 1/^2), where the various elements refer to the amplitudes in the separate minima, while the 
wave functions correspond to (l/-\/2, —1, l/-\/2, 0) and (0, 1/^/2,-1, l/-\/2). These latter are higher in energy 
by a factor 

= = ,3.) 

As 5 0, Ai? vanishes, leading to a static, but disordered JT phase, with the octahedra equally likely to 
have any of the four tilt distortions. The resulting lattice symmetry is pseudo-cubic, in that all orientations 
are equally likely. This four well approximation becomes exact as — s- 00, and Figure 3 shows that it remains 
qualitatively correct for all values of v. The solid lines in Fig. 3 are the exact eigenvalues of Eq. 32a, as 
discussed in Appendix V. The axes of Fig. 3 are in normalized units, introduced in Appendix V: i/ = ^^/q, 
E± = Aaa. 

There appears to be an asymmetry in the above wavefunctions, since they are all centered on wells 2 or 
3. However, linear combinations of these wavefunctions can be generated which are centered on the other 
wells. For example, \/2(0, l/\/2, 1, l/\^) - {l/s/2, 1, l/\/2,0) = (-l/\^,0, l/s/2, 1). 

7b. Phase Transitions: Dynamic JT Phases 

As the temperature is lowered, intercell coupling will lead to an ordered low temperature phase, with 

< cos(f> > or < siruf) > non zero. A complete solution of Eq. 29 is prohibitively difficult, particularly since 
most of the parameters are not well known. In this section, a number of simplifications are introduced to 
make the problem more tractable, allowing a determination of the conditions under which an LTO (dynamic 
JT) phase might be stable. These simplifications are: (1) since a key question will be to establish whether 
an LTO-like phase can be stabilized in the absence of strong umklapp scattering, it will be assumed that 
7*^' = 0. (2) It will be assumed that < cos(f) > and < sincj) > are small near the transition (i.e., that the 
transition is second order, or weakly first order), so that H2 can be linearized in these quantities. 
(3) Even with these assumptions, remains complicated: 



[r^ - ro+ + 2ro- + sFoo] d4^_„3^ ^ ^ , „-^3^ 

+ [de^ - {Too + -^)R'^]R'^{cos(f> < coscj) > +sin(j) < sin(f> >) 

+Ci^'^P-E?'{cos(f> < coscj) > —sincj) < sincf) >). (38a) 
In the same spirit of replacing Eq. 30 by Eq. 31, this may be replaced by the simpler form 

H2 = Ho{cos(j) < coscj) > +sin(j) < siruj) >), (386) 
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which involves only one unknown parameter, Hq. To see how good this approximation is, consider the case 
in which the tilting Hamiltonian, Huh = + ff| of Eq. 29, is dominated by the vHs JT effect (i.e. the 
terms involving af^). In this case, the electronic term can be integrated out as in Eq. 30, leading to 

Huh = 2 — |(cos0+ < cos(f> >)^ — {sin(f)+ < siruf) >)^|. (39) 

Figure 4 shows that even for this highly singular potential, the simpler form of Eq. 31 plus 38b provides 
a reasonable approximation. In approximating Eq. 39, the parameter Hq must be chosen to be negative; 
hence, only this regime will be explored in detail below. 

With the intercell coupling given by Eq. 38b, the perturbation matrix may be recalculated. The 
Hamiltonian matrix, Eq. 35, becomes 



/ x + AHi y + {A-B)H2 -{y+{A + B)H2) 

y + {B-A)H2 x + BHi y+{A + B)H2 

y-{A + B)H2 X - AHi y + {B- A)H2 

\-{y-{A + B)H2) y + {A-B)H2 x-BHi 



(39) 



where x = Hn -E, y = H12 - SE, A =< coscj) >, B =< sin(p >, Hi = HQexp{-l/8i^), and H2 = HiS/V^. 
For arbitrary < sincf) >^ 0, Eq. 39 must be diagonalized numerically. However, for the LTT phase (< 
sin(j) >= 0), the solutions can be foimd analytically. The potential wells arc labelled in such a way that the 
< cos(j) > terms lower the state (1,0,0,0) (state tpi) and raise the state (0,0, 1,0) (f/'a) in energy. Then, in 
the LTT phase {B = 0), two solutions have ipi = 0, with energies 



2SH - AHi/2 ± J2H^ + A^{H^/4 - 2ff|(l - 2S^)) - 2SHAHi , ^ 

-^1,2 = ^ _ ; V^Oa) 



while the other two have V's = and 



2SH + AH1I2 ± ^2H-^ + A^{Hll^ - 2HI{1 - 2S^)) + 2SHAHi 
EsA = . (406) 

In the above equations, E is measured from Hn, and H = HnS — H12. Figure 5 illustrates how the energies 
in the four wells vary as a function of < coscj) >, for several values of S, both for the LTT phase (< sin(j> >= 
- dashed lines) and for the LTO phase (< siruf) >=< coscf) > solid lines). [Note that there is a break in 
slope in the LTT phase at (< cos(j) >)^ = 1/2, since the constraint (< sin(j) >)^ + (< coscf) >)^ < 1 comes 
into play. In the figure, it is assumed that < sincf) >= ^1 — (< coscfi in this case. However, this regime 
has no immediate relevance, since < cos<f) > is always found to be < 0.7 in the numerical calculations for 
the LTO phase.] 

Prom Figure 5, the energy is always lowered when < coscp >^ 0, and at low enough temperatures 
there will be a transition to a dynamic JT phase. In mean field theory, the phase diagram may readily be 
calculated from the self-consistency condition for < cos(f> >: 

Trcoscbe-"/^'''^ 

< '^"i' >= Tre-H/^^BT ' (41) 

where Tr stands for the trace. In the four-state model, coscj) vanishes on average in wells 2 and 4, and 
coscf) = ±6-^/^" in wells and 3(-). Thus 

< coscf, >- EM.-P^.y-'^^'^^ (42) 

^coi,cp^ e ^-Ei/ksT p-E2/kBT f^-Es/kBT p-Ei/ksT' V"*^^ 

where Pij is the probability that well i is occupied in the state of energy Ej . Note that any term common to 
all the Ei^s cancels out of the ratio in Eq. 42. Thus, Eq. 42 depends on three parameters, S (or, equivalently 
u), H, and Hq. At sufficiently high temperatures, Eq. 42 has no non-zero solutions. As T is lowered, there is 
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a critical temperature below which there are values of < coscf) >^ 0. Figure 6 illustrates the solutiou for the 
LTT phase, plotting the right-hand side of Eq. 42 against < coscp > at a number of different temperatures. 
The solutions of Eq. 42 correspond to the intersection of the solid and dashed lines. These solutions lead to 
the evolution of the ordered phases, < coscf) > (T) shown in Figure 7. 

7c. Phase Transitions: LTO vs LTT Phase 

Figure 7 compares the resulting phase diagrams for both the LTT and the LTO phases, for several sets 
of parameters. The parameters were chosen to approximately match the critical temperature of the LTO 
phase in La2Cu04. Figure 8 shows the free energies of the two phases, 

for the same sets of parameters. 

From the above figures, it is possible to understand the competition between the LTT and LTO phases 
as a competition between energy and entropy. First, consider the small overlap case. Fig. 5a. As 5 — > 0, the 
LTO phase becomes two independent LTT-type transitions, in wells 1 and 2, so two branches of the energy 
curve are lowered by virtually the same amount as the single branch in the LTT phase. Thus, for a fixed 
value of < cos(f) >, the LTO phase has a greater entropy associated with it. However, this greater entropy 
in turn means that at a fixed temperature, the self-consistent value of < coscf) > will be smaller in the LTO 
phase than in the LTT phase, Eq. 41. Hence, for small S values, the LTT phase is always energetically 
preferred. 

As S increases, however, the factor of level repulsion adds a further stabilization to the LTO phase 
(Fig 5b-d). The splitting of the originally degenerate pair of LTO levels drives one lower in energy than 
the corresponding LTT level. Also, the mixing of different wells brings two LTT levels fairly close together, 
so the entropy differences between the LTO and LTT phases are reduced. The result is that both phases 
are close in free energy over most of the temperature range. Fig. 8, and that both phases allow non-zero 
solutions for < cos<f) > out to significantly higher temperatures as S increases. However, the LTO phase wins 
out at the higher temperatures, and over a wider temperature range as S gets larger - the phase diagram 
is illustrated in Fig. 9. Other choices for the parameters H and Hq would mainly shift the scales of the 
phase diagram, without altering its fundamental character. Note that the LTT phase never increases its 
transition temperature significantly above the S limit, whereas the LTO phase turns on at higher T's 
with increasing S. 

While the above behavior is qualitatively what might be expected, it should be cautioned that the 
above calculations are perturbative in 5", so the detailed nature of the large S results should be treated with 
caution (for instance, the LTO transition temperature appears to diverge as 1/2)- A more detailed 

calculation would require the inclusion of more distant overlaps in Eq. 39. However, the reduction of the 
Mathieu problem to a four-state problem remains valid, as long as 4q: is large compared to /jbT. Moreover, 
as Fig. 9 illustrates, for some parameter choices the LTO phase remains stable down to small S'-values. 

From Fig. 7, the transition from the disordered JT phase (HTT) to the LTT phase is second order, 
with < coscf) > acting as an order parameter, but the LTO^HTT and LTO^LTT phase transitions are first 
order, with discontinuous jmups in < coscp >. Note that in the LTO phase, for S > 0.1, < coscf) >~ 0.4 
at low temperatures - i.e. , the macroscopic average tilt is noticeably smaller than its instantaneous value, 
-\/ < coscf) >^ -I- < sincf) >2 ~ 0.56. 

While the present dynamic JT calculation has been carried out for the high-Tc cuprates, it is interesting 
to compare it with calculations of the multiple phase transitions in the ferroelectric perovskite BaTi03[38]. 
There, a similar competition between energy and entropy was found, with the largest energy lowering as- 
sociated with the static (rhombohedral) JT phase, and increasingly large entropy contributions in the case 
of two, four, or eight potential minima being involved in the dynamic JT effect. However, in the simplest 
calculation, all four phases had a transition at the same temperature, so the static JT phase was stable at all 
temperatures. By assuming that the interaction energy was different in the different phases, it was possible 
to reproduce the observed sequence of phase transitions. In the present calculation, the overlap parameter 
naturally provides the stabilization for the dynamic phases, with only the static phase stable as S ^ 0. 

8. Discussion 
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8a. Dynamic JT 

Figure 9 is the chief result of this paper. It confirms the suggestion[10,ll] that the LTO phase can be 
described as a dynamic JT phase involving a splitting of the vHs degeneracy. 

The present paper has contrasted two possible origins of the LTO phase, as a static JT phase stabilized 
by umklapp scattering (Section 5), and as a dynamic JT phase, Section 7c, which is naturally related to the 
LTT phase. By comparing the gap functions of the two models, Figs. 2 vs 7a, it can be seen that, for small 
S values, the dynamic JT effect can mimic the effect of a static transition driven by umklapp scattering, 
thereby confirming the speculation made in Paper VIIIA[10]. For larger S values, the situation is more 
complicated, with the LTO transition becoming first order. 

Much work remains to be done, particularly in simultaneously accounting for both the electronic and 
structural aspects of the transition, and in understanding how electronic properties are modified within the 
dynamic JT phase (antifcrromagnetism, superconductivity, ...). Before this can be done, however, some 
more fundamental questions must be answered, such as, what does the Fermi surface (or even the Brillouin 
zone) m,ean in a dynamic JT system, where the local symmetry is not the same as the global symmetry, and 
indeed where the local symmetry can fluctuate in space and time. 

8b. SDW-CDW Competition Revisited 

This paper has explored possible generalizations of CDW's and SDW's in the presence of a vHs. While 
the simple CDW can describe the LTT phase, it seems to be ruled out by a symmetry argument in the LTO 
phase (Ref. [9] and Appendix III). Nevertheless, there arc at least two mechanisms by which vHs-related 
effects could stabilize an LTO phase. Which of these mechanisms is actually operative in LSCO is a question 
which requires considerably more research. Nevertheless, it seems appropriate to point out that even at the 
molecular level, there is often a competition between JT and spin-orbit cffects[12,13], and this competition 
is reminiscent of the usual CDW-spin-density wave (SDW) competition in lower-dimensional metals in the 
presence of a peak in the dos. 

Near half filling, correlation effects become important in destabilizing CDW-like phases. On the other 
hand, there has been clear experimental[39] and theoretical [40] evidence for the importance of spin-orbit 
coupling in the undoped materials, and specifically in the Neel antiferromagnetic or spiral magnetic phases. 
This can be understood, in the context of the present paper, as follows. To split the vHs degeneracy requires 
breaking the degeneracy of the two Cu's in the orthorhombic unit cell. This can be done if they have opposite 
spins, but this requires spin-orbit coupling to modify the electronic bands. 

I earlier suggested that the phase diagram of LSCO resembled a crossover from SDW-like behavior near 
half filling (the antiferromagnetic state) to 'incipient CDW-like' behavior, including Peierls distortion, as 
hole doping is increased [28]. The present results suggest the following modification: the hole-doped LTO 
phase is stabilized by the dynamic JT effect, whereas near half filling spin-orbit effects are more important 
(due to on-site Coulomb repulsion), and may stabilize a static LTO phase. This phase may be related to 
the proposed 'flux phase' [18]. Whereas in principle, a single spin-orbit coupled phase could persist for all 
dopings, there is considerable experimental evidence for a transition between two phases as a function of 
doping, with a (nanoscale) two-phase regime between the vHs and half filling [2 8]. 

In Section 2, it was briefly noted that c-axis dispersion acted in the same way as a CDW to split the 
vHs degeneracy. It would be interesting to study in more detail whether changes in intcrlayer hopping could 
drive a phase transition, particularly in light of Anderson's ideas that interlayer coupling plays a special role 
in stabilizing high-Tc superconductivity [41]. 

9. 'Excitons' vs 'Cageons' 

When I initiated this series of calculations on the vlls[42], a major premise was that inter- vHs scattering 
could promote large electron-phonon coupling via vHs nesting [N.B. not conventional nesting], leading to 
short-range CDW order. Indeed, the term 'exciton' was introduced to point out that the strong scattering 
is associated with electron-like sections of one vHs scattering off of hole-like sections of the second vHs. The 
CDW order in this case would be an 'excitonic instability' analogous to the spin-density wave instability of 
chromium. 

The 'excitonic' properties of the model were discussed in paper IV[43]; the CDW in V[27]. The present 
series of papers, VIIIA-C, are an extension of V: there is a structural instability, but it is not simply 
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dcscribablc as a CDW (the formalism of V is still relevant for extending the present model to a quasi-two- 
dimensional system and incorporating mode-mode coupling). 

Under these circumstances, the term 'excitonic' does not seem to be particularly suitable for describing 
the system, since the strong electron-phonon coupling will lead to excitations closer to polarons - or even 
solitons, as discussed below. Hence, I am introducing the term 'cageon', which is intended to better describe 
the JT excitations of electrons coupled to tilts of the octahedral CuOg cages. 

Going beyond the mean field calculation, the corner sharing of the CuOg octahedra suggests that there 
should be long chains of LTT phase, and that defects must be introduced into the chains. A plausible model 
would be to have islands of LTT phase separated by LTO-like domain walls, which would switch an x-directed 
domain into a y-directed one. Such domain walls would behave as solitons, and the doping dependence of 
LBCO could be interpreted in terms of the generation of these solitons. Thus, at x = 0.12 (6% of the La 
replaced by Ba), the material is in a pure LTT phase. As x is reduced, the material transforms to the LTO 
phase. This transformation could be accomplished via soliton generation - as x is decreased, the density of 
LTO solitons increases. The octahedral shear ei2 would simply be proportional to the soliton density. This 
idea will be pursued further in a future publication. 

In particular, it has been found that the degeneracy of the vHs's is split in the solitonic model. It is 
just this splitting which underlies the 'valence bond density wave' calculations of Ref. [10]. Hence, it seems 
likely that the results of that paper hold for the dynamic JT model of the LTO phase - in particular, the 
phase diagram of the transitions HTT — > LTO — > superconductor should continue to hold in the dynamic 
JT model. 

I would like to thank J. Zak for explaining the role of the Berry phase in the JT effect. Publication 545 
from the Barnett Institute. 

Appendix I: Renormalized Microscopic JT Band Structure 

The present analysis, in terms of p±. is very convenient, but how can it be related to a microscopic 
Hamiltonian, such as that developed in VIIIB, Appendix I? In particular, the analysis of Section 5 requires 
being able to define a\{l), a\{l) - i.e., on each atomic site. By contrast, the microscopic Hamiltonian involves 
interatomic hopping between Cu and O atoms. In this Appendix, I show that a site model can be derived 
from the microscopic Hamiltonian, if the sites are not individual atoms, but clusters of atoms. 

The simplest cluster is a single octahedron - or more simply, a Cu atom with the four surrounding 
planar O's, since the present model does not incorporate either the Cu dz'^ nor the apical O p orbitals. 
However, one octahedron is too small. There is only a single Cu-O antibonding level (twofold degenerate 
due to spin) per Cu atom, and hence no JT degeneracy. Hence, the appropriate cluster contains a square of 
four octahedra. The antibonding 'band' contains four levels, the middle two of which are degenerate, in the 
absence of strain or tilt coupling. This is a JT degeneracy, since the antibonding 'band' is half filled (the 
model of VIIIB contains no 0-0 hopping, so the vHs falls exactly at half filling). 

The dispersion of the four levels can readily be recovered from VIIIB. 

E=^± ^Ji^r + ^W, ilia) 

W = tl^^^[cos'e^sin'i^) + 0lsinH^cos\^)] + tl^^^[cosXsin\^) + ^IsinHyCos^^)]. (716) 

Here, is the splitting of the Cu and O levels, and tcuOx (tcuOy) and 9x [Oy) are the Cu-O hopping energy 
and octahedral tilt along the x (y) axis, respectively. For a 'sample' two cells by two cells, kxa/2 and kya/2 
are restricted to the values and it/2, leading to four possible values of E, Eq. Ila, with eigenfunctions 
shown in Fig. 10 (top). (Strictly speaking, these are not small clusters, because periodic boundary conditions 
are assumed. Calculations on real c;lusters would reveal eigenstates of the same overall symmetry as in Fig. 
10, but with modified eigenvalues.) 

This same procedure may readily be generalized to larger 'samples'. For a sample N x N cells, there are 
N"^ levels in the antibonding band, of which N states are degenerate (when the strains and tilts are absent) 
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at the vHs energy level - those states for which kx + ky = ir/a. For these states, the energy can be written 
as Eq. II, with 

W = Wa + WbCos{k^a), (/2a) 



e+ + 



4(/?^-l) 




3ma^ 



4(/3^ + l) 



(Rl-Rl)] 



(72c) 



3ma^ 



where r* , P^r, Ri = 3ma^sin?9i/8, and t^^MOi = ^Cuoi^ ~ Sae^j/r*), as discussed in VIIIB (see also the last 
paragraph of this appendix). This equation offers a convenient microscopic estimate of the electron-phonon 
coupling constants. 

In order to find the renormalized bands corresponding to pij of Eq. 4, it is convenient to look at 
how larger clusters are built up. Fig. 10 (bottom) shows the eigenfunctions associated with the 4x4 
clusters. Beneath each figure, the corresponding eigenvalue is indicated in the form ij. For simplicity, 
this notation corresponds to the untilted clusters {dx = 0y — 0), for which W in Eq. lib is given by 
4:W = if^cuOx j^CuOy "^^^ 4x4 eigenfunctions are built up from the 2 solutions by repeated tiling of the 
4x4 cell by the 2x2 cell or its negative. Since there is no intercell mixing, the Cu-0 band separates into a 
superposition of four overlapping but noninteracting subbands. This pattern holds for larger cells, but with 
some minor complications: (a) all of the '+' Cu atoms need not have the same amplitude [example: for the 

one-dimensional chain 8 cells long, the eigenfunction (+ + + H ) is really (a, b, b, a, —a, —b, —b, —a), 

with a/b = \f2 — 1]; and (b) most levels are degenerate in pairs, leading to eigenvalues of more complicated 
form. These complications do not affect the subband separation. 

From studying larger cells, a separate energy dispersion can be determined for each subband. These 
energy dispersions have the form of Eq.Il, with a restricted range of k^^ ky, for each subband. Thus, for 
subband I (IV), kx and ky must both be less than (greater than) 7r/2, while for subband II (III), kx (ky) is 
greater, while ky {kx) is less than tt/2. This restriction can most elegantly be carried out by introducing new 
artificial Brillouin zone boundaries, as illustrated by the dashed lines in the inset to Fig. 11. For the case of 
untilted molecules, the subband dispersions are illustrated in Fig. 11. It can be seen that subbands II and 
III are degenerate (in the untilted case), and overlap the vHs, while band I(IV) lies below (resp. above) the 
vHs. For the purposes of the present calculation, bands I and IV may be neglected, while bands II and III 
correspond to the electronic bands discussed in the text, e.g., Eq. 4. In particular, since these bands are 
asymmetric in kx vs ky, their degeneracy will be lifted by a nonvanishing LTT-type tilt. 

It should be noted that Schulz[17] has introduced a similar formalism with (spin-dependent) creation 
operators associated with each vHs. 

ADDENDUM to Paper VIIIB [11]. In Paper VIIIB, the covalent overlap enhancement factor 
was introduced as 



Here and r^r are related to the Slater-Koster[44] parameters of p-d overlap of cr or tt symmetry as follows: 
Ta = ypda, T-w = —VpiiTr (N.B., Vpdn is negative). At the time of writing, I was unable to find theoretical 
values for the two overlaps separately. Now Grant and McMahan[45] have provided ab initio calculations 
for tetragonal La2Cu04 of sufficient detail to allow an estimate of these parameters. Their calculations 
provide two estimates of r^^. From McMahan and Grant, Table I: t{dx2-y2,p^) = y/^Ta/^ = lA7eV; 
t{d^22_^2,p„) = Ta = 0.50ey. These estimates are not exactly equal, since the wave functions are Wannier 
functions which are not pure p and d states, due to overlap with higher orbitals. From Grant, Table 2.5 
t{dxy,P%) = = 0.72eV. Averaging the estimates, this leads to a theoretical value for ~ 1.6, 
somewhat smaller than the value used in Ref. [11]. This value is in good agreement with recent cluster 
calculations [46]: Vpda = l-5eV, Vpd-n = —0.7eV, yielding = 1.5. I would like to thank A. McMahan for 
providing me with copies of Ref. [45] . 



Appendix II: Phonon-mediated Electron-electron Interaction 
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The strain terms can be decoupled from the Hamiltonian, Eq. 1, by defining 

e±(9)=e±(g) + ^, {Ilia) 

ei2 (g) = ei2 (g) + ^ . (7/16) 

The reduced Hamiltonian becomes H = H'^^^ + H'_^ + H'_ + Hq + Hph + H^, with 

H'str = ^(c+e+(-g)e+(g) +C_e_(-g)e_((?) +2C66gi2(-g)e~i2(g)), (//2a) 

9 

H'o = EK^(|)(y^o'^(t) + rvi2(-g))) - E YPi2(-g)pi2W, (//2c) 

9 9 

where r[,± = -/^|/2C±, T'^ = -FiJlGCee, = al-G±F±/2C±, Y' = Y -GmFee/^Cee, J'± = Gi/2C±, 
and Jq = G%q/2C%q. The term H'^^^ is decoupled from the remaining terms, and since it is purely quadratic, 
does not contain any interesting dynamics, hence can be neglected. 

The form of //pft,, Eq. 5d, does not allow a similar decoupling of the tilt terms. Such a separation can 
be obtained for a related 'toy' Hamiltonian, which consists of Eq. 1, with Hph replaced by 

9 

If fo = (To + rg)/16, fa = Ta/ie, and oj'^' = Cj'^{qR), then Eq. 113 will be identical to Eq. 5d for the 
soft mode, q = 0, (f = qn = (7r/a, 0) or (0, 7r/a). Equation HS can now be incorporated into Eqs. II2b,c 
by substituting Tq^. Tq^ = Fq + Tq^. and r'2 = r2 + F2 into the latter equations, and replacing 

Hph with 

H'p^ = u;''nf{o). {im) 

At this point, the tilt couplings can be formally decoupled from the toy Hamiltonian by a procedure similar 
to the strain decoupling, by defining 

mi)=mi)+^^, {lua) 



urn 



In this case, H = H'^^^ + Hr + Hp, with 

9 "i=± 

'j+ I N , , J- , , , , Jo 



//h =a;^'^f (0) + E(E[ro.^i(T)^^(|)J + T^o^T^^^^IO' 

9 i=± 

= - E(t''+(~^^''+(«^ + ^P-(-g)p-(g) + ^Pi2(-?)pi2(g)), (//56) 



with 



r<2 ^e/2 
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and 



/066 J- 2 



Even for the toy model, this separation must be treated with caution, because of the unusual form of the 

substitution, Eq. 114. In the strain case, the separation can be carried out by a canonical transformation, 
whereas in the present problem, this does not seem to be the case, due to the quadratic term in R in Eq. 
114. 



Appendix III: Structural Instabilities in the Uniform LTO Phase 



This Appendix analyzes the question of whether a uniform LTO phase could be brought about by 
electron-phonon interaction. That is, can the large dos associated with the vHs's be driven below the Fermi 
level by a lattice distortion. It has been shown that the LTO transition does not split the degeneracy of 
the vHs's[9], but in itself, this is not sufficient. There are two alternative means by which the dos could be 
shifted below the Fermi level, without splitting the vHs degeneracy. Here it is shown that neither of these 
effects arises in LSCO, in the absence of spin-orbit interaction. 

These effects are: (1) there is a single peak in the dos, due to both (degenerate) vHs's, but the structural 
transition shifts this peak below the Fermi level; or (2) umklapp scattering couples two vHs's opening a gap 
in the dos, as discussed in Section 5. Case (1), which will be discussed in Section Ilia, would arise if, for 
example, due to the orthorhombic distortion, the vHs's no longer occured at the corners of the orthorhombic 
Brillouin zone, or if the transition reduced the magnitude of the average 0-0 hopping matrix element. 



Ilia: Orthorhombic Distortion 



The Brillouin zone of Fig. lb for the LTO phase has been oversimplified by the neglect of the orthorhom- 
bic splitting. The LTO phase modifies the Brillouin zone of the HTT tetragonal phase in two ways. First, 
the principal axes a* and b* are rotated by 45° with respect to the tetragonal axes, a and b, Fig. lb. (The 
real space cell is doubled in area, so the Brillouin zone is halved.) Secondly, there is a small orthorhombic 
distortion, a* 7^ 6*, which was neglected in Fig. lb. 

This distortion can readily be incorporated into the tight-binding calculations of Appendix I. When 
these calculations are repeated for the larger LTO unit cell, two changes arise. First, each band becomes 
two- fold degenerate (due to the two Cu's per unit cell). Secondly, in all dispersion relations, such as Eqs. 
II, 12, the following substitutions must be made: 




{Illla) 



kya^^-^^, {I I lib) 

with i5+ = A;+a*, i5_ = k-b*, and k± the wave vectors along the new principal axes. As a result of this 
transformation, the topology of the Fermi surfaces is maintained as the Brillouin zone is stretched into a 
rectangle. In particular, the vHs's continue to intersect the Brillouin zone boundaries in the corners of the 
zone, as in Fig. lb. Thus, the orthorhombic distortion does not alter the ratio of the area of a given Fermi 
surface to the total Brillouin zone area - and hence does not change the doping at which the vHs's coincide 
with the Fermi level. 

This still does not rule out the possibility of a shift of the vHs away from the Fermi level. Thus, the 
doping X = Xc at which the vHs coincides with the Fermi surface is controlled by the curvature of the 
Fermi surface at the vHs. The ratio of the area of the hole Fermi surface to the total Brillouin zone area 
is (1 -I- Xc)/2. But the curvature of the Fermi surface is proportional to the 0-0 hopping parameter, too- 
when too = 0, the Fermi surface is square. Since the orthorhombic distortion changes all 0-0 distances 
(and produces two inequivalent distances), it could lead to a change in the average value of too, and hence a 
shift of Xc- Such an effect is likely to be small, since one 0-0 separation increases while the other decreases, 
so the corresponding changes in too tend to cancel. 
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Illb: Umklapp Scattering and Group Theory 



The role of umklapp scattering can best be appreciated by analyzing a Fermi surface away from the 
vHs, on the side of overdoping. Fig. 12a. Simple zone folding from the tetragonal to the orthorhombic 
Brillouin zone produces the lens-shaped orbits shown in Fig. 12a. However, these orbits are anomalous, 
having discontinuous slopes at the Brillouin zone boundary. In most band structures, umklapp scattering 
causes an interaction between carriers from opposite sides of the Brillouin zone, opening a gap between 
successive bands. This caiises the Fermi surfaces of each band to intersect the Brillouin zone boundaries 
orthogonally, so that in crossing the zone boundary (in an extended-zone scheme) the carrier stays within 
the same band, and there is no slope discontinuity. 

Such an interaction will arise in a tight-binding model in the following fashion. The orthorhombic unit 
cell contains twice as many atoms as the original tetragonal cell. This will double the size of the energy 
eigenvalue matrix, and the number of bands, leaving the bands degenerate in pairs. To split the degeneracy 
of the bands requires that the matrix elements of formerly equivalent atoms (e.g., the two Cu atoms in the 
cell) be different. However, there is a (glide) symmetry operation which can translate one Cu into the other, 
so the matrix elements can only differ in overall sign ~ e.g., terms proportional to ±sin6, in Appendix I. 
Since these terms enter the eigenvalue equation only in the square, they cannot produce an energy splitting. 
Thus, umklapp coupling is absent - pi2=0. 

But how can the slope discontinuities of the lens orbits be accounted for? This can best be understood 
by a general symmetry argument, based on the Bmab space group of the LTO phase. There is a glide 
plane perpendicular to the b* direction, which interchages the two Cu atoms in a cell. Because of this, all 
eigenfunctions on the a*-face of the Brillouin zone must be doubly degenerate[47]. In the presence of such 
degeneracy, the slope of the E{k) curves can be discontinuous. Such a situation always arises in the presence 
of glide planes, and is perhaps best known for the hexagonal close packed space group[48]. In this case, one 
zone boundary does not introduce a gap, and the E{k) curve from the first band merges continuously into 
the second band on crossing the zone boundary. It is often convenient to ignore the zone boundary, and 
work with a larger k-space zone[49]. Such a double zone for the LTO phase is illustrated in Fig. 12b: the 
zone is doubled in the X* direction (along a*), while a gap is allowed along the Y* direction. Thus, lens 
orbits appear only along the zone face at Y*, with open orbits along X*[50]. 

IIIc: Spin-Orbit Interaction 

Spin-orbit coupling lifts most of these degeneracies, so that the enlarged k-space zone can no longer be 
used[51]. The Bmab space group is equivalent to the Cmca space group, which includes the structure of 
the elements Br, I, and Ga. The group theoretical analysis for Ga has been carried out by Koster[52], who 
finds that spin-orbit coupling eliminates most of the degeneracy on the faces of the Brillouin zone. This is 
illustrated in Fig. 13. A complication arises in that the orthorhombic unit cell, with axes a*, b* , and c is not 
primitive. Figure 13a shows this cell, with the primitive cell (containing only half as many atoms) inscribed 
in it. The Figure also shows the Brillouin zones corresponding to the primitive cell (Fig. 13b) and the 
orthorhombic cell (Fig. 13c) - for convenience, the latter zone will be called the pseudozone. Its importance 
arises because it is the natural zone for tight-binding model calculations, particularly in the two-dimensional 
limit when the energy bands are assumed independent of kz. Figure 13 illustrates how the true zone may 
be folded into the pseudozone. Figure 13b shows the first pseudozone inscribed in the full zone, while Fig. 
13d shows how the second pseudozone is reassembled from the leftover parts of the full zone. This folding 
produces two bands, which are degenerate within a tightbinding calculation such as that of Appendix I. Note 
the relative orientation of the two pseudozones in the full zone: the second zone is predominantly displaced 
from the first by the Q-vector Z* = Z/2. This is because the unit cell (Fig. 13a) includes contributions 
from two Cu02 planes displaced along the c-axis, which are equivalent in the tight-binding scheme. 

By including additional terms in the tight-binding calculation, it is in principle possible to couple the 
two bands and remove the twofold degeneracy. However, a group theoretical analysis shows that not all of 
the degeneracy can be lifted. In Figs. 13b, c, the hatched areas show the regions in the Brillouin zone (all 
confined to the surface of the zone) in which the wave functions are two-fold degenerate (neglecting spin 
degeneracy) in the absence of spin-orbit coupling: this happens throughout the L — X — N plane and along 



22 



the line M N. When spin-orbit coupUng is included, the degeneracy is lifted except on the regions which 
are heavily shaded: along the X — L line and at the M point. 

Prom the pseudozone of Fig. 13c, the approximate two-dimensional Brillouin zone (Figs, lb, 12a) is 
found by neglecting the band dispersion along the c-axis, which should be a good approximation for the 
cuprates. From Figure 12, it can be seen that spin-orbit coupling lifts the degeneracy over must of the 
Brillouin zone boundary. In particular, the lens orbits will be continuous at the zone boundary, due to the 
opening of a gap, and spin-orbit coupling can lead to an umklapp gap even at the vHs, thereby stabilizing 
a uniform LTO phase, as discussed in Section 5. 

Appendix IV. Pseudospin Formalism 

The calculations of Sections 6b, c are similar to previous studies of structural phase transitions[12-14,29], 
which have introduced pseudospin formalisms, to describe either the electronic or the phonon modes. These 
formalisms arc often based on an imperfect analogy, but are valuable because spin systems are understood 
so much better than nonlinear phonon systems. To clarify comparison with these works, the formalism is 
discussed in this Appendix. The pseudospin formalism applies to the order-disorder limit of structural phase 
transitions. That is, it is assumed that at each lattice site the system is disordered into one of two (or more) 
possible configurations. In the present problem, an electronic pseudospin can be introduced, with the 'up' 
spin associated with occupancy of the pn vHs, and the 'down' spin with p22 electrons. Then 

a^i = al^aii - al^a2i, {IVla) 

(^xi = a\-a2i + 4.aii, {IVlb) 
ayi = iial^au - 0^02,), {IVlc) 
where the subscript i labels the atomic site. The spin operators satisfy the commutation relations 

The p operators can then be rewritten in terms of pseudospins. In particular, 

/^-(O) = ~ 4fc«2;s) = '^(^zi = Sz, {IV3a) 

k i 
i 

Thus, condensation into the LTT phase corresponds to a pseudospin ferromagnetic transition within the 
plane, < Sz >^ 0; the overall tetragonal symmetry arises because the interlayer coupling is antiferromagnetic. 

In the present model, the pseudospins are again only of approximate validity. The problem lies in the 
term J2ik ^kdlkf^ik- In an order-disorder model, each electronic state would have to be either in state 1, 
associated with the X-point vHs, or state 2. associated with the Y-point vHs, in which case a,\i.aik + a\i^a2k = 
1 for each fc-value. In fact, while the dos peaks are split, there remains a significant overlap of the two bands, 
and for these states, the sum can reach a value 2. Hence, the Eh-teim must be retained in Eq. 19a, leading 
to a very different form of gap equation from a pseudospin calculation. 

Appendix V: Solutions of Mathieu's Equation (Eq. 32a) 

The solutions of Eq. 32a can be written in the form of a Fourier series[31] 



e= J2 ''rne'"'', {VI) 



m= — oo 



where the boundary conditions require that m -I- 1/2 is an integer. Substituting VI into Eq. 32a yields the 
recursion relation ^ 

ami^m - am^) + -P{am-4 + am+i) = 0, (V2) 
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where £m is the eigenvalue and a = fi^ /2R^ . Eq. V2 only mixes m-values separated by ±4, so there are four 
families of solution, depending on whether the series VI contains terms with m equal to 3/2, 1/2, -1/2, or 
-3/2. Since the recursion relation contains only rri^ , these terms are degenerate in pairs for any value of (3. 
When /? = 0, each term in VI is an exact solution, with eigenvalue = am^ . Nonzero (3 leads to mixing 
of these states, and for large /3, each level is fourfold degenerate (corresponding to independent oscillations 
about one of the four potential minima). 

In finding the eigenvalues of Eq. 32a, it is convenient to generalize the equation to: 

+ (a' - 2q'cos{lz'))y = 0, (F3) 

where the number of minima, I, is arbitrary, and the boundary condition is that y{z' + I-k) — —y{z'). This 
agrees with Eq. 32a when a' = En/a, q' = (3/2, z' — (f), and 1 = 4. In turn, Eq. V3 can be reduced to 
the canonical form of a Mathieu function of fractional order[33] by the substitution z = lz'/2, a = Aa'/P, 
q = M'/l^, so 

-4 + (a - 2qco.s{2z))y = 0. (V4) 

From the boundary condition, the solution can be written in the form of series VI, with m = n+p/l, where 
n is an arbitrary integer and p is an odd integer < I. These solutions are the fractional Mathieu functions 
of order p/l, and the eigenvalues can be read off of Fig. 11 of McLachlan[33] (p. 98). For fixed n, there are 
I solutions lying between the solutions of the integral Mathieu functions, of orders n and n + 1 (actually, in 
the 'stable' zone between the solutions a„ and 6„+i) - see Fig. 14 for I = 3,4,5. For even I, the solutions 
are all doubly degenerate, while for odd I, there is an additional singly degenerate level whenever p = I. As 
I increases, the allowed states fill the 'stable' intervals of Mathieu's equation, with gaps between successive 
n values, much as the band structure of a one-dimensional metal fills in as more and more atoms are added 
to the chain. 

Specializing now to the case of interest, I = 4, the eigenvalues are shown in Fig. 14, constructed by 
interpolation from McLachlan's Fig. 11. The nature of the eigenstates can be appreciated by going to the 
(7 = limit. In this case, the eigenstates are ypc = cos{p6/2) or yps = sin{p6/2), withp = 1, 3. Thus, y^^ = 1, 
0.5, 0, or 0.5 in the wells 1, 2, 3, and 4, for either value of p. The difference is that y^c has additional maxima 
outside of the potential minima. The functions yps are similar, with wells 1 and 3 interchanged. Since yic and 
yis are degenerate, a number of alternative solutions can be constructed, including ones centered on wells 2 
or 4. As q increases, the interwell tunneling probability decreases, and for sufficiently large barriers, the wave 
functions should localize within a single well. It might be anticipated that the individual wave functions 
would narrow, causing the overlap of wave functions between wells to decrease. A plausible measure of this 
decrease would be i^p = yp^(2) / yp^{l) - i.e., the overlap probability of the wave function being found at 
the center of well 2 given that its peak value is centered in a neighboring well, 1. Surprisingly, r]p = 0.5, 
independent of q. Localization arises from the mixing of the p = 3 states into the p = 1 states, so that the 
gap AE between these two levels can be taken as a delocalization energy. 

The eigenvalues may be found more precisely (Fig. 3) by deriving a continued-fraction eigenvalue 
equation. Letting Vnp = am/am-4, with m = {p + 8n)/2, Eq. V2 can be rewritten 

q 



a-{p + 8n)2 - qvn+i,p ' 



iV5) 



Here p is a positive integer which labels the four series, p = 1,3,5,7. Equation V5 is readily solved as a 
continued fraction 

. (^6) 



a - [p + 8n]^ — 

a-[p + 8(n+l)f 



a - [p + 8(n + 2)]2 - ... 

Equations V2 and consequently V6 hold for n > 1, but the equation for n = is more complicated, since it 
mixes the series for p and p = 8 — p. This equation can be written in the form 

a-q-y-p"^ 

vip = , {V7a) 
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with 7 = apjija^ji. Equating the right-hand side of Eq. 6, for n = 1, to the right-hand side of Eq. 7 gives 
the eigenvalue equation for the Mathieu function - actually a pair of equations for v\p and v\f. The quantity 
7 can be found by subtracting the two equations. Letting 

_ -y^ ^ qjvip - Vip) 
y~ 2q 

then 

7± = -2/ ± a/1 + 2/^- 

The two possible 7 values give the two degenerate eigenstates for each energy. Substituting either into Eq. 
V7a yields a single eigenvalue equation for a{q), which is solved numerically to generate Fig. 3. 
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Figure Captions 

Fig. 1 Fermi surfaces at vHs associated with (a) LTT phase; (b) LTO phase. 

Fig. 2 LTO phase diagram based on strong umklapp scattering (Eq. 22). Solid lines correspond to Jo = 200, 
300, 400, or 500K, in order of increasing gap, D. 

Fig. 3 Lowest energy eigenvalues of Mathieu's equation, Eq. 32a. Both the exact solutions (solid lines) and 
the approximate solutions of Eq. 36d (£'+ = dashed line; E- = dotted line) are shown. 

Fig. 4 Sohd hnes = potential well of Eq. 39, E = 2Htiit/Q.''lB? , with < cos(f> > = 0.2, for (a) LTT and (b) 
LTO phases. Dashed lines are approximations to the potential using Eqs. 31 + 38b. 

Fig. 5 Eigenfunctions of Eq. 39 for the LTT (dashed lines) and LTO phases (solid lines), assuming H = 400K, 
Hly = -200K, for several values of the overlap S = 0.001 (a), 0.01 (b), 0.1 (c), and 0.5 (d). 

Fig. 6 Graphical solution of Eq. 42, where Arg is the right hand side of Eq. 42 (solid lines) and the dashed 

line is < cost/) >. The solution of Eq. 42 is at the point where solid and dashed lines intersect (filled 
circles). The calculations assume the same parameters as Fig. 5, with S = 0.01, and the various solid 
lines correspond to T = 100, 80, 60, 40, or 20K, in order of increasing < cos(f) >. 

Fig. 7 Equilibrium values of < cos(j) > (solutions of Eq. 42) for the LTT (dashed lines) and LTO phases (solid 
lines). The parameters of Eq. 42 were chosen as H = 400K, Hq = -200K, with variable S. In order of 
decreasing values of < coscf) > or Tq, the values of S in part (a) are 0.0001, 0.01, 0.02, 0.05, and 0.1. In 
part (b), the values of S are, in order of increasing Tq, 0.3, 0.4, 0.48, 0.5, 0.52, and 0.55. The transition 
temperature Tc has a minimum around 5 ~ 0.2. 

Fig. 8 Free energies of LTT (dashed lines) and LTO phases (solid lines), for the same parameters as in Fig. 7. 
Increasing S corresponds to increasing magnitude of F in part (a), but decreasing magnitude of F in 
part (b). Filled circles in part (b) show LTO-LTT crossover. 

Fig. 9 Phase diagram of dynamic JT transitions. Solid lines show HTT— >LTO and LTO— >LTT phase bound- 
aries, Tc{S), for the same parameters as Figs. 5-8: H = AOOK, H'q = Dashed lines are for the 
alternative parameters H = 600if , H'q — ~9Q0K . 

Fig. 10 Eigenfunctions of Eq. II, for 2x2 (top) and 4x4 (bottom) clusters. In each cluster, only the Cu's are 
indicated. The ±-sign corresponds to the relative phase of the Cu dx2_y2 orbitals (periodic boundary 
conditions assumed). 

Fig. 11 Partial densities-of-states for the four subbands corresponding to Fig. 10a (dashed lines) along with the 
total dos (solid lines). Inset: Brillouin zone showing subband boundaries. 

Fig. 12 Fermi surfaces of LTO phase, with orthorhombic distortion: (a) ordinary cell, including spin-orbit 
coupling. Circles on zone boundary show points where bands must be degenerate, filled circle implies 
degeneracy for all points along kz, open circle implies degeneracy for special points only, (b) Double 
zone appropriate when spin-orbit coupling is absent. 

Fig. 13 Group theoretical analysis of orthorhombic Brillouin zone, (a) Real space cell of orthorhombic phase, 
showing orthorhombic cell (long dashed lines) and inscribed primitive cell (short dashed lines). Dotted 
lines = glide planes. The open and filled circles represent Cu atoms, with filled circles representing 
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atoms at a; = 0, a*, and open circles atoms at x = a*/2. The solid lines represent the canted planes 
of the 'planar' O's. (b) Full Brillouin zone, corresponding to primitive cell. Dashed lines = inscribed 
pseudocell. Light shading represents portions of Brillouin zone surface where the bands are twofold 
degenerate in the absence of spin-orbit coupling (i.e., the X — L — N face and the line — M). Heavy 
shading indicates those portions where the degeneracy persists in the presence of spin-orbit coupling 
(the line L — X and the point M). (c) Orthorhombic pseudozone. Shading has same meaning as in 
Fig. 13b. (d) Zone folding construction of second pseudozone. Here shading illustrates how sections are 
folded from the full zone into the pseudozone. 



Fig. 14 Eigenvalues of fractional Mathieu's function. Solid lines = solutions to the integral Mathicu's function 
and boundaries of stability for the fractional values. Fractional solutions are shown for orders 1/3 (short 
dashed lines), 1/4 (long dashed lines), and 1/5 (dotted lines). [After Ref. [33].] 
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